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Abstract 

We show that the nearest-neighbor spacing distribution for a model that consists of random 
points uniformly distributed on a self-similar fractal is the Brody distribution of random matrix 
theory. In the usual context of Hamiltonian systems, the Brody parameter does not have a definite 
physical meaning, but in the model considered here, the Brody parameter is actually the fractal 
dimension. Exploiting this result, we introduce a new model for a crossover transition between 
Poisson and Wigner statistics: random points on a continuous family of self-similar curves with 
fractal dimension between 1 and 2. The implications to quantum chaos are discussed, and a 
connection to conservative classical chaos is introduced. 
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It is well known that the spectral statistics of almost all quantum systems, whose classical 
analogs are chaotic, are described quantitatively by Gaussian ensembles of random matrices 



flflQQ. 



Perhaps less well known is the fact that the eigenvalue spacing distribution for 
the Gaussian orthogonal ensemble of 2 x 2 random matrices also describes homogeneous 
Poisson point processes in R 2 |3j. The nearest-neighbor spacing distribution (NNSD) of 
random points uniformly distributed on a line is given by 

P P (S) = exp(-S), (1) 

and the NNSD of random points uniformly distributed on a plane is given by [3] 

Pw{S) = \s ex P (-^ 2 ), (2) 

which is, in fact, the Wigner distribution of random matrix theory (RMT) 5]. (In Eqs. (JIJ 
and (J2J), S is a dimensionless scaled spacing.) Evidently, random points on a line are uncorre- 
lated, whereas random points on a plane are actually correlated (in the sense that the points 
tend to avoid each other). Even so, the latter result is reasonable since there is an additional 
degree of freedom that allows the points to spread out. Based on this intuitive interdepen- 
dence between "point repulsion" and dimensionality, and in strict analogy to energy-level 
statistics, it is tempting to indiscriminately conjecture that the nearest-neighbor statistics of 
random points on a fractal set with noninteger dimension between 1 and 2 are described by 
an intermediate distribution in-between Poisson and Wigner. In fact, this conjecture turns 
out to be correct as we shall demonstrate below. 

Our interest in the nearest-neighbor statistics of random points on fractals was sparked 
by the provocative results ((TJ) and (J2J), and the prospect that random points on fractals could 
be conceptualized as new models of intermediate statistics (between Poisson and Wigner). 
We however are not the first authors to consider the statistical properties of fractal sets. 
Two decades ago, Badii and Politi 6f studied (for completely different reasons) the nearest- 
neighbor distance distribution of random points on a strange attractor. These authors also 
obtained the probability distribution of nearest-neighbor distances 5 among N points chosen 
randomly and uniformly on a Cantor set with capacity dimension D . Interestingly, they 
noted (quite tersely) that the asymptotic distribution (appropriate for large N) J(| 

P BP {5, N) = 2D N(25) D °- 1 exp [-N(26) D °] , (3) 
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could be "recognized as a Brody distribution" . Clearly, Eq. (J3J) is not the Brody distribution 
0, but we show that when the nearest-neighbor distance is rescaled by the mean nearest- 
neighbor distance (which is a contrivance familiar to practitioners of RMT) the result is 
indeed the Brody distribution. More generally speaking, we show that the Brody distribu- 
tion is the NNSD for points selected uniformly at random from a self-similar set K C R d 
with similarity dimension d s > 1, and that the Brody parameter is the (relative) similarity 
dimension of K (i.e. d s — 1). (In this paper, d is the Euclidean space dimension.) The goals 
of this paper are to derive this result, to introduce a new model for a Poisson-to-Wigner 
crossover transition based on this result, and to explicate the physical implications. The 
derivation is simple and direct; the result itself is far more interesting than its proof. A 
discussion of the physical implications is deferred to the conclusion. 

We begin first with the derivation of the spacing distribution. Suppose that iV points of 
a self-similar set K C H d are chosen randomly and uniformly. The probability P(s)ds of 
finding the nearest neighbor to a given point at a distance between s and s + ds is equal to 
the probability of finding one of the (N — 1) points at a distance between s and s + ds from 
the given point and the (N — 2) remaining points at a distance greater than s. Let &{s) 
denote the probability of finding a point within a distance s of a given point. The probability 
of finding one point at a distance greater than s is then (1 — &(s)), and for (N — 2) points, 
the probabilities are multiplicative due to the fact that all points are chosen independently. 
Thus, 

P(s)ds = 7 1 (l-$(s)) {N - 2) dp(s) } (4) 
accounts for the fact that the nearest neighbor could be 



where the prefactor rj = (N — 1 
any one of the (N — 1) points [2jJ, and da>(s) = <s\s)ds is the probability of finding a 
point in a shell with inner and outer radii s and s + ds centered about the given point. 
The probability of finding multiple nearest neighbors is ignored since the probability of such 
an event is higher order in ds and is therefore insignificant compared to the probability of 
finding a single nearest neighbor. 

It now remains to specify the probability 2>(s). Clearly, the typical number of neighbors 
of a given point will vary more rapidly with distance from that point as the dimension 
increases. The probability (p(s) is, by definition, the ratio of the number of points within 
some prescribed distance to the total number of points, that is, 2>(s) = N(s)/N(Jl), where 1Z 
is the radius of the rf-dimensional ball that contains all iV points. The number function N(s) 



can be directly obtained from the so-called "mass-radius scaling law" for fractals (see page 
40 of Mandelbrot's book @): M(r) = M(R){r /TZ) Dm = \M(K) /K Dm ]r Dm . In this formula, 
M(r) and M(1Z) are the masses contained within balls of radii r and TZ, respectively, and 
D m is the mass dimension. For regular (i.e. strictly self-similar) fractals, D m — d s . It might 
seem peculiar to speak of masses here, but it is equivalent to the concept of numbers of 
points within balls of a specified radius if an individual sample point becomes the unit of 
mass. We can thus equitably think of the above mass law as a "number-radius scaling law" . 
So, the probability of finding a point within a distance s of a given point is governed by the 
power law 

<e{s) = As d % (5) 

where the coefficient A = l/lZ ds , and d s > 1 is the similarity dimension of K. (Note that the 
probability t(s) is a unitless number since A has units of 1/ (length) ds .) Therefore, Eq. (J1J) 
becomes 

P{s)ds = n (1 - As d °) {N ~ 2) Ad s s da - l ds. (6) 

Recall that dtp(s) = & f (s)ds = Ad s s ds ~ 1 ds. It is straightforward to verify that the proba- 
bility density P(s) is already normalized (i.e. P(s)ds = 1). We could now consider the 
arg e iV limit of Eq. (JHJ), and in doing so, we can invoke the so-called Poisson approximation 



22|, and this gives the asymptotic probability density 

P(s) = NAdsS^' 1 exp (-NAs da ) as iV -»■ oo. (7) 

Equation (J7J) is essentially the distribution obtained by Badii and Politi in 1985 (c.f. Eq. (J3J) 
and note that Dq in their formula is equivalent to d s in Eq. ( |7| l). 
Next, we calculate the mean spacing s = J™ sP(s)ds: 



pi 

7Z(N - 1)4 / u d ° (1 - u d °) {N ~ 2) du 
Jo 

K(N-l) [ v l ^(l-v^ N -^dv 



Jo 

= n{N-l)B{\ + 1/4, N- 1) 

= ^r(iv)r (i + 1/4) /r (n + 1/4) . 

In the second line, we have made a change of variables u = s/lZ, and in the third line, 
we have made one further change of variables v = u da . The integral in the third line we 



recognize as the Beta function B(fi, v) with parameters fi = 1 + l/d s and v = N — 1, and 
this then gives the fourth line. We then used the usual relation between the Gamma and 
Beta functions to arrive at the fifth line. It can be shown that the term 

m) - 1 L 0fiY) asiV^eo, (8) 



T(N+l/d s ) N 1 /** V 
and therefore, the asymptotic mean nearest-neighbor spacing is, to leading order, 

T — - — as N — > oo. (9) 



N 1 /^ V d s 

Introducing the rescaled spacing S = s/s and taking the limit N — > oo, the distribution 
P(s) in Eq. (0) becomes the distribution 

P B (S; q = d s -l) = adsS^" 1 exp (-aS ds ) , (10a) 

where 

(Notice that the rescaling of s by s was doubly beneficial; both explicit dependences on 
1Z and iV in Eq. (J7J) have been removed.) The distribution Pr(S\ q) [Eq. (fTUj)] is, in fact, 
the Brody distribution [7] with Brody parameter q equal to d s — 1 [3]. We refer to the 
number d s — 1 as the relative similarity dimension of K since this number is the difference 
between the similarity dimension of K and the similarity dimension of a line (the simplest 
self-similar object) which is equal to 1. Equation (jlOJ) is valid for random points on any 
self-similar subset of H d with similarity dimension d s > 1. If K is a classical (nonfractal) 
self-similar set (i.e. a rf-dimensional cube), then d s = d and Eq. (fTU|) reduces to the NNSD 
for a homogeneous Poisson point process in H d (see Ref. |3(). 

In studies of quantum chaos, the Brody distribution has sometimes been used (as a purely 
phenomenological distribution) to describe the nearest-neighbor energy-level statistics of 
quantum systems that undergo a direct transition from Poisson-like to Wigner-like statistics 
as a system parameter is varied. (A classic example is the diamagnetic Kepler system Jj|.) 
In the present context, a Poisson-to-Wigner transition can be realized by considering point 
processes on a family of self-similar sets whose dimension ranges between 1 and 2 as some 
set parameter is varied. In actual fact, we are introducing a new model (that does not 
involve random matrices) for a Poisson-to-Wigner crossover transition, and this model is 
special since the intermediate statistics are described exactly by the Brody distribution. As 



a concrete example, we now study point processes on the family of Koch fractals in R 2 . 
These fractals can be thought of as the attractors of a one-parameter family of iterated 
function systems (IFSs). The similarity transformations defining the IFS involve a rotation 
which is conveniently parametrized by the angle 9. When 9 = 0, the attractor is a line, 
and when 9 = tt/2, the attractor is the famous Sierpinski-Knopp plane-filling curve, whose 
image is a solid isosceles triangle in R 2 1(3] . For intermediate values [i.e. 9 G (0, 7r/2)], the 
attractors are various self-similar curves of prescribed dimension d s e (1,2). The nearest- 
neighbor statistics of the random points undergo a continuous transition from Poisson to 
Wigner statistics (see Fig. [TJ as the self-similar set continuously deforms from a line to a 
plane-filling curve (i.e. as the rotation angle 9 varies between and 7r/2). 

For clarity, we mention some pertinent numerical details. Random points on these 
fractals were selected using the random iteration algorithm (RIA) 11]. The distance be- 
tween a given point Xj and its nearest neighbor is defined (using the Euclidean metric) by 
si = min ly/ixi - Xj) 2 + {yi - yj) 2 ^ for i, j = 1,...,N (j ^ i). Although {sj define a 
set of spacings, the NNSD is actually defined in terms of the scaled spacings Si — Si/s, 
where s = 4 J2iLi * s the (numerically-calculated) mean nearest-neighbor spacing. We 
constructed the histograms by binning all values of Si and then normalizing the area under 
the histogram to unity. Each histogram was constructed from one sample of N = 20 000 
random points. We mention here that by using all of the points selected by the RIA in 
the statistical analysis, we introduce some error due to finite-size or edge effects. These 
errors are statistically insignificant as long as iV is sufficiently large (see Eq. (J0|)). Although 
not absolutely necessary, we use the Levenberg-Marquardt method to determine the 
numerical value of the parameter d s that gives the optimal fit of the Brody distribution 
Pb{S;cI s — 1) to the histograms. This number df can then be immediately compared to 
the theoretical value. The purpose of this procedure is to test how accurately the Brody 
distribution [Eq. (jlOj) ] reproduces the histogram data obtained from particular realizations 
of the model. Of course, each realization (in general) yields a unique histogram (and hence 
a unique df), and so it is more informative to average over several (say n) realizations and 
subsequently define df = df ± a, where df is the average df value obtained from the n 
realizations and a is the standard deviation. The percentage error (denoted by e) of df 
(obtained from n = 10 independent realizations) relative to d s is typically under 1%. We 
have in fact studied point processes on many of the well-known classical fractals in R 2 , R 3 , 
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FIG. 1: A crossover transition between Poisson and Wigner statistics resulting from point processes 
on the family of Koch fractals in R 2 . The left panel shows random points on several fractals that 
belong to the family (each one specified by a particular choice of the rotation angle 6). The 
exact similarity dimension d s of each of the fractals (given to three decimal places) is indicated 
on each window. The right panel shows the corresponding NNSD of the points on each fractal. 
The numerical data (i.e. df = df ± a and e) indicated on each window are as described in the 
text. The dotted, dashed, and solid curves are the distributions Pp(S), Pw(S), and Pb{S; df — 1), 
respectively. 
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and R 4 [24], and in almost all cases the accuracy is comparable. 

The nearest-neighbor energy-level statistics of some quantum systems execute a Poisson- 
to-Wigner transition as the underlying classical dynamics monotonically change from being 
completely integrable to completely chaotic. One example is the family of Robnik billiards 
Q|. A monotonic transition between integrabililty and chaos is however quite exceptional. 
More typically, the degree of "chaoticity" of the classical dynamics (as measured by the 
volume fraction of phase space filled with chaotic trajectories) changes in a complicated way 
as a system parameter is varied monotonically, and so the energy-level statistics will not un- 
dergo a direct transition from Poisson to Wigner. For example, the energy-level statistics of 
the hydrogen atom in a van der Waals potential undergo a Wigner-Poisson-Brody-Poisson- 
Brody-Poisson- Wigner transition as the appropriate system parameter is monotonically var- 
ied in a specified range Regardless, in the intermediate regime between integrability 
and hard chaos, the Brody distribution (albeit a pure surmisal) has often been found to 
be a good delineation of the energy-level spacing histogram. There are, in fact, systems 
for which the statistical confidence is high (an example is the ripple billiard This is 

not to say the Brody distribution is now established as a distribution that quantitatively 
describes energy-level statistics in the intermediate regime, but rather that after 30 years of 
pervasive use with considerable success, the Brody distribution has garnered an undeniable 
phenomenological significance. (Of course, other distributions have been proposed and used 
to interpolate between the Poisson and Wigner limits; we cite here a few of the more popu- 
lar distributions |3, (ill (3 till Q • These efforts cannot be disregarded, but they are not 
directly relevant to the present discussion.) Given our present result and the phenomeno- 
logical status of the Brody distribution in studies of quantum chaos, there is a profound 
implication that transpires: Phenomenologically , the energy levels of a typical time-reversal 
invariant quantum system, whose dynamics in the classical limit are mixed, have the same 
nearest-neighbor statistics as random points on a fractal with dimension in-between 1 and 
2. This phenomenological corollary, together with analogous corollaries in the limiting cases 
of integrability [Eq. (£[])] and hard chaos [Eq. (J2*|) ]. offer a new phenomenology for quantum 
chaos: the nearest-neighbor energy-level statistics of a typical time-reversal invariant quan- 
tum Hamiltonian follow the statistics of (i) random points on a (one-dimensional) line if the 
classical limit is integrable; (ii) random points on a (two-dimensional) plane if the classical 
limit is fully chaotic; and (iii) random points on a fractal set with dimension in-between 1 
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and 2 if the classical limit is mixed. 

This phenomenological behavior is quite puzzling. Why should the energy levels of a 
quantum Hamiltonian behave (insofar as their nearest-neighbor statistics) like random points 
on a fractal or on a plane? This is a very difficult question to answer since there is no direct 
connection between point processes and quantum mechanics. Point-process models (PPMs) 
and random-matrix models are both stochastic models, but unlike random-matrix models, 
PPMs do not inherently contain any of the structure of quantum mechanics, and so it 
is difficult to understand why point-process statistics should have any relation to energy- 
level statistics. Surreptitiously, the fundamental link is classical mechanics. The model 
of random points on a fractal can be regarded as a simple stochastic model for chaotic 
dynamics on a Poincare section. In mixed Hamiltonian systems, regular and chaotic regions 
are comingled, and the chaotic regions, in particular, are fractal in nature. If we restrict 
our scope (at least initially) to two-degree-of-freedom billiard systems (such as the family 
of Robnik billiards), then we know that the chaotic trajectories explore (in a seemingly 
random fashion) a fractal subset of the Poincare section having dimension in-between 1 and 
2. Clearly, the NNSD of the "chaotic points" on the section (corresponding to a chaotic 
trajectory) must be Poisson-like in the near-integrable regime, Wigner-like in the chaotic 
regime, and some intermediate distribution in-between Poisson and Wigner in the mixed 
regime. The intermediate distribution must also have built in point repulsion. If random 
points on fractal sets embedded in R 2 really are apt models of Hamiltonian chaos (on 
a Poincare section), then the intermediate distribution should be the Brody distribution. 
If so, then there is an even deeper corollary: Phenomenologically, the energy levels of a 
typical time-reversal invariant quantum system (whose classical analog is nonintegrable) have 
the same nearest-neighbor statistics as the chaotic trajectories of the underlying classical 
Hamiltonian. Of course, we have not explicitly demonstrated that PPMs correctly describe 
the nearest-neighbor statistics of chaotic trajectories, and we can only begin to do so through 
numerical experiments. This shall be the subject of a future paper. The purpose of the above 
discussion was merely to introduce the idea of linking PPMs with classical mechanics, and 
to outline one of the potential implications. For the present, we must settle for the less 
fundamental, but nonetheless far-reaching precursor (italicized above). 
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